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FOREWORD 


This is a report on the research topic “Aerodynamic Drag Reduction.” This work 
was conducted as part of the research activities of the Institute for Computational and 
Applied Mechanics (ICAM) which is now a subprogram of the Institute for Scientific and 
Educational Technology (ISET). Within the guidelines of the research project, special 
attention was directed to “Numerical Investigation of Flow Around Rectangular Cylinders 
with and without Jets.” The period of performance of the specific research ended December 

31, 1998. 
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NUMERICAL INVESTIGATION OF FLOW AROUND 
RECTANGULAR CYLINDERS WITH AND WITHOUT JETS 

By 

S. N. Tiwari 1 and S. B. Pidugu 2 
Department of Mechanical Engineering 
Old Dominion University, Norfolk, Virginia-23529 

SUMMARY 

The problem of flow past bluff bodies was studied extensively in the past. The problem of 
drag reduction is very important in many high speed flow applications. Considerable work has 
been done in this subject area in case of circular cylinders. The present study attempts to investi- 
gate the feasibility of drag reduction on a rectangular cylinder by flow injection by flow injection 
from the rear stagnation region. The physical problem is modeled as two-dimensional body and 
numerical analysis is carried out with and without trailing jets. A commercial code is used for this 
purpose. Unsteady computation is performed in case of rectangular cylinders with no trailing jets 
where as steady state computation is performed when jet is introduced. It is found that drag can be 

reduced by introducing jets with small intensity in rear stagnation region of the rectangular cylin- 
ders. 
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Chapter 1 

INTRODUCTION 


The problem of flow past bluff bodies is one of the oldest in the fluid mechanics, and it 
remains one the most important problem due to the complexity of the fluid dynamic phenomena 
despite their geometrical simplicity. Cylinders of different cross sections such as circular, rectan- 
gular and their combinations have been used in the past, due to many practical applications. Vari- 
ous fluid dynamic issues of prime interest and importance in this problem are drag minimization, 
vortex shedding and aeroelastic instability. 

The subject of bluff body flows has recently been receiving a great deal of attention. This 
is largely because of its importance for energy conservation. For example, road vehicles must now 
meet more stringent fuel consumption requirements which translate into the need for reduced 
drag. Large structures such as tall buildings must also be designed so as to minimize convective 
heat losses. Besides, large structures must be designed so as to avoid potentially disastrous wind 
induced large amplitude oscillations. This needs knowledge of the vortex shedding characteristics 
and is also crucial in the design of vortex shedding flowmeter for measuring flow inside closed 
conduits. Therefore, this field of study is highly desirable. 

Roshko[l] extensively conducted wind tunnel experiments to investigate the drag and 
shedding frequency of two dimensional bluff bodies. A correlation between bluff bodies of differ- 
ent shapes was obtained in order to calculate the drag by combining free streamline theory with 
experimental results. Bostock and Mair [2] experimentally determined drag and base pressure 

coefficients around rectangular and D-shaped cylinders for Reynolds number between 0.5 X 10 5 
to 1.6 X 10 5 . For cylinders of D-shaped section, it was shown that the pressure distribution on the 
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curved surface and drag are considerably affected by the state of boundary layer at the separation 
point, as for a circular cylinder. Furthermore, they showed empirical correlations for drag based 
on the known values for the constituent half-bodies are, in general, unsatisfactory. Bearman and 
Truman [3] measured base pressure coefficient, drag coefficient and Strohaul number of rectangu- 
lar cylinders for Reynolds number between 2 X 10 4 to 7 X 10 4 . It was also demonstrated that the 
flow around sections is found to be strongly influenced by the percentage of the trailing edge cor- 
ners. Knauess and John [4] obtained a correlation between Reynolds number and Strauhal number 
for elliptical cylinders through experiments. A correlation of data on Strouhal number variation 
with angle of incidence, based on the projected height of the model, is obtained for both square 
and elliptical models in the range of Reynolds numbers between 300 and 1200. Strouhal numbers 
of rectangular cylinders have been determined experimentally by Atsushi [5] as a function of Rey- 
nolds number (70 - 2 X 10 4 ) and aspect ratio(l to 4). They observed, for a certain range of Rey- 
nolds number for cylinder with aspect ratio of 2 and 3, flow pattern changes abruptly with a 
sudden discontinuity in the Strouhal number. 

Various numerical investigations have been also conducted on the flow past bluff bodies[6- 
8]. Clements [6] developed an inviscid model of two dimensional vortex shedding in order to 
numerically predict flow features behind bluff bodies, form of vortex shedding and Strouhal num- 
ber. Davis and Moore [7] presented numerical solutions for two dimensional flow about rectan- 
gles in infinite domains. The range of Reynolds numbers used was from 100 to 2800. The 
numerical work was extended to include the effects of confining walls [8]. 

Many investigators used different shapes of bluff bodies as well as different methods of 
flow control in order to make their studies useful for practical applications. Jet blowing method 
near the rear stagnation region is widely used in order to reduce the drag by various investigators 
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[9,10], A splitter plate behind the cylinder is a popular technique used by many investigators to 
reduce the drag [11]. Both methods have advantages and disadvantages. The advantage of using a 
splitter plates is that no power is required to reduce drag significantly. This can be done easily by 
placing a plate in the existing system without need of any complicated mechanisms. However, this 
method is not suitable for cases where rectangular cylinders have non-zero angle of attack. Any 
mechanism to make splitter plate to adjust its direction in the flow brings additional mechanical 
complications, thus making the system lose its simplicity. Moreover, it is difficult to control the 
amount of the effect of the plate. On the other hand, jet blowing method is capable to change 
direction as well as its strength; and, more importantly, the jet can be actuated only when it is 
needed without any geometrical change. Particularly, this method is very useful when compressed 
air is readily available. The present study attempts to investigate the feasibility of drag reduction 
on a rectangular cylinder by flow injection from the rear stagnation region. 

As mentioned earlier, many investigations are limited to flow past circular and rectangular 
cylinders. Considerable amount of work also has been done by introducing jets as well as splitter 
plates m the rear stagnation point, in case of circular cylinders. But there is little or no information 
available on the effect of jet injection from the rear stagnation point of rectangular cylinders on 
the flow characteristics. The present study is directed to the investigation of the effect of free 
stream Reynolds number and aspect ratio on the flow past rectangular cylinders with and without 
jets. Heat transfer aspects are not considered in this study. 
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Chapter 2 


GENERAL FORMULATION 


2.1 Physical Model 

The physical model, as shown in Fig. 2.1, consists of a rectangular cylinder placed half- 
way between two channels separated by a distance L. The rectangular shape has known length A 
and height B. The diameter of the orifice on the rear face of the rectangle is kept at 10 percent of 
the height of the rectangle when jet injection is considered. The channel length is also sufficiently 
long in order to exclude the wall end effects that might arise in the computation. The upstream 
face of the rectangle is located at a distance x r from the inlet. The problem is limited to two 
dimensions in order to reduce the computational time and memory requirements. Free stream 
Reynolds number based on rectangle length is varied from 1 X 10 4 to 1.25 X 10^. 


2.2 Numerical Formulation 

2.2.1 Governing equations 

The two-dimensional time dependent continuity and momentum equations for a Newto- 
nian fluid with constant density and viscosity are given by 
Continuity 


9u 9v 
9x + 9y 


= 0 


X-momentum 


( 2 . 1 ) 


l; ( p u)+ s (pu2, + | (puv) = _ ^ + ^ 


A 


^ 9 2 u 9 ju 

v 9x 2 9y 2 J 


( 2 . 2 ) 
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Fig. 2. 1 Physical model for two dimensional flow past a rectangle cylinder with jet 


5 



Y-momentum 


d , . 0 , . 0 , 2 . 0p 

dt P V + 5x * >UV + 5y P V ~ - dy + P 


f- j 2 ->2 \ 

0 V 0 V 

2 + ri 


0x 0y' 


y 


(2.3) 


Equations (2.1)- (2.3) can be written in a compact form by denoting the dependent variable as <|> . 
Thus, the governing differential equations can be expressed in general form as 


^(p<t>) + V(pv<(>) = V[r (|> (grad<t))] + S (t , 


(2.4) 


where the corresponding values for <{>,r ^ , and are listed in Table 2.1. 

2.2.2 Boundary conditions 

Boundary conditions required for the physical problem are shown in Fig. 2.2. No slip con- 
ditions are imposed on the confining walls as well as on the surface of the cylinder. Neuman con- 
ditions are employed at the downstream outlet where gradient of the dependent variable in the 
normal direction is fixed at zero. Free stream velocity is given at the inlet. When jet is introduced, 
the fixed velocity is given at the jet slot. 

2.2.3 Solution Procedure 

Detailed formulations and solution procedures are given in Fluent manual [12,13]. Here, 
few necessary details on the theory are provided for the analysis of the present problem. The gov- 
erning equations, Eqs. (2.4), are reduced to finite difference analogs by integration over the com- 
putational domain dividing it into number of cells. 
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Table 2.1 Values of i n Eq. (2.4) 

Equation <J, T + 

mass 

x-momentum 


y-momentum 











u = 0 
v = 0 


Outlet 


'///////////////////////////////////////////////yzzzzzzz//zzzzzzzzL^^. 



3u 

3x 


dv 

3x 


3p 

3x 
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Fig. 2.2 Boundary conditions 



u = 0 
v = 0 
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Time integration of the transient terms can be descritized as 


t + At 

J (p < t>)dt = (p^r 1 - (p^j 1 

t 

where n is the time level 


(2.5) 


Time integration of other terms of Eqs. (2.4) can also be written, in general form, as 

t + At 

J (pu<(>)dt = [f (pu<>f + 1 +(1 -fj (pu^fjAt 


t 


(2.6) 


where f is a factor which determines whether the solution is explicit or implicit. Fluent assumes 
=1, which corresponds to the case of fully implicit; hence 

t + At 


J (pu<(>)dt = f (pu^j 1 + 1 At 


(2.7) 


f 


t 

Divergence theorem is employed to perform volume integration of the space divided over 
the control volume as 


(P v< t ) ) dv - /(pV^ldA (2.8) 

V A 

where V is the volume of the control volume and A is the area of the control surface. After descri- 
bing the integral on the control volume as shown in Fig. 2.3, the surface derivative terms in gov- 
erning equations become (for example in the x-direction) 
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(pu<j>A) c - ((pu<>A) w + (p u4> A) n - (pu<()A) s ) 



, (4>„-4>p) 

0 Ax 



+ S a A 


where source term can be expressed as = -(p e - p w ) 


(2.9) 


For two dimensional case, we can write complete descritized analog of governing equation 
as 


ip«v)r'-(P4>y>: 

At 


>n- 1 


+ (pu<(>)“ 1 - (pu<t>)” + 1 + (pu<t>A)" + 1 - (pu<|>A) 


n + 1 


. n + 1 


= n 


r r (4’ e -'fp) i" + l r (4> p -» w ) -|" +l 

L * A Je _ L r * St A J W 


L ♦ Ay A J n L 4_ Ay' A J S A (210) 

By collecting terms with respect to space, Eq. (2. 10) can be written in algebraic form as 


Vp = M °' d <l>p ld + a e <t> e + a w <)> w + a n <}) n + a s <|> s + S, 


( 2 . 11 ) 


where M old = (pV) old which corresponds to n-level in time in Eq. (2.1 1), while other terms are 
evaluated at (n+1) level in time. Finally, Eq. (2.1 1) can be written in compact form as 
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(2.12) 


V»p = S^nb + VM^V 

nb 

where subscript nb refers to neighboring points. 


The equations are solved using the SIMPLE algorithm which is based on using a relation- 
ship between velocity and pressure correction in order to recast the continuity equation in terms ot 
pressure correction calculations. 

Using in Eq. (2.12), the momentum equations are expressed as 


x-direction: a p u p = ^a nb u nb + (p w - p e )A + ( M°' d <!) p ld ) x 

nb 

(2.13) 

y-direction: a p v p = ^a nb v nb + (p n - p s ) A + ( ^f ld (() p ld ) y 

nb 

(2.14) 


In general, velocities obtained from Eqs. (2.13) and (2.14) do not satisfy the continuity 
equation until pressure values are substituted in the source term of the equation. To overcome this 
problem, the continuity equation is used to obtain the correct pressure values. Thus, in turn, cor- 
rect velocities are obtained which satisfy the continuity equation,. This is the essence of the SIM- 
PLE algorithm. 

* 

The SIMPLE algorithm starts with the substitution of a guessed pressure field, p , into 

momentum equations which can be solved to obtain guessed velocity field, u . These guessed val- 
ues are related as 


U P = U P + U P 


V = V* + V 

p p p 


Pe = Pe +P C 


(2.15) 

(2.16) 
(2.17) 
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where prime denotes corrected values. Substitutions of Eqs. (2.15)-(2.17) into Eqs. (2.13) and 
(2.14) with corrected values, yield 

x -direction : a p u p ' = £a nb u' nb +(p w -p e )A (2.18) 

nb 


y-direction: a p v p ' = J>nb v nb +(p n -p s )A (2.19) 

nb 

Equations (2.13) and (2.14) are used to relate the pressure and velocity corrections as 
x-direction: u p ' = J-(p' w - p ' e ) A (2.20) 

p 


y-direction: v' = — (p' n -p' s )A 


(2.21) 


In Eqs. (2.20) and (2.21), the terms containing the neighbor influences are dropped for conve- 
nience. However, these terms vanish when convergence is achieved. 

To obtain correction terms of pressure, the descritized continuity equation is employed 
which can be expressed as 


(puA) e -(puA) w + (pvA) n -(pvA) s = 0 


( 2 . 22 ) 


In Eq. (2.22), u and v are replaced by starred (*) and primed (‘) values using Eqs. 
(2.15) and (2.16). Moreover, the primed values are represented by the pressure correction terms 
through Eqs. (2.20) and (2.21). Finally the pressure correction equation can be expressed as 

(pu*A) e - (pu*A) w + (pv*A) n - (pv*A) s 


+ (PA) e/A (Pp Pp ) (P A )w>'a \ (Pw'P'd) 


(Ap ) e -P -E 


W ( A P>w 


+ (pA)n (X^ (p ; -P , N)-(P A )s ( -^y-(P / s-P / p)= 0 


(2.23) 
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The Eq. (2.23) is used for the pressure correction which is then used to compute velocity correc- 
tion via. Eqs. (2.20) and (2.21). Finally, the velocity and pressure are updated using Eqs. (2.15)- 
(2. 17). The system of implicit equations are solved by using the line Gauss Siedal method. 

2.2.4 Time Step 

Physical problems which are time dependent should be solved by choosing proper time 
step in order to capture physical phenomena accurately (time accurate). Since, the Fluent code 
uses fully implicit scheme, and there is no way to find out the CFL number condition in the Fluent 
code, several numerical experiments have been conducted to check the effect of the time step on 
the computations. After several trials, different values of time steps are selected for different inlet 
velocity conditions. 

When the jet is introduced in the wake region, the problem is solved by using the steady 
state assumption. This assumption is mainly based on the justifications given in the work done by 
Atsuchi [13]. Steady state assumption reduces the computational time (by several times) as well 
as large memory requirements to store the data at different time steps. 

2.2.5 Grid System 

Several grid sizes have been tested to compare the results and performances. A grid size of 
1 15 X 62 is employed for computational purposes as shown in Fig. 2.4. To get the physical phe- 
nomena correctly near the surface region, grid points are highly clustered near the wall. When the 
jet is introduced in the wake region, a more dense grid of size 130 X 71 is employed. In this case, 
grid is also clustered near the jet mixing region as shown in Fig. 2.5. 
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Fig. 2.4a Grid system around a rectangular cylinder without jet 



Fig. 2.4b Grid system around a rectangular cylinder without jet 
(in the vicinity of the wall) 
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Fig. 2.5a Grid system around a rectangular cylinder with jet 



Fig. 2.5b Grid system around a rectangular cylinder with jet 
(in the vicinity of the wall) 
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2.3 Important Parameters 

Important parameters are presented in this section. 

2.3.1 Base pressure coefficient, Cp 



P-P 

U 

2 r °° 


00 


2 

00 


2.3.2 Strouhal number, St 

St 

U oo 

where f is the frequency of vortex shedding. The Strouhal number is a very important non-dimen- 
sional number for the flow around bluff bodies which shed vortices into the down stream. The 
variation in the flow pattern can be related to the change in the value of the Strouhal number. 

2.3.3 Drag coefficient, Cd 

C = D 
d 1 2 

ip U ooh 

2 00 

2.3.4 Lift coefficient, Cl 

C. = — 

1 1 2 

2P00U ood 

Here D and L are drag and lift forces. 


2.3.5 Jet blowing momentum coefficient, Cm 

To evaluate and normalize the strength of the jet, the jet blowing momentum coefficient, 
C m is used as given in [10], i.e.. 
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c 


m 


Pi d i u ? 

1 2 u 
2P~ u ooB 


In this study, Pj and are same and, therefore, the strength of the jet is normalized with 

respect to the cylinder height and free stream velocity. Here subscripts °° and j represent free 
stream and jet conditions respectively. 
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Chapter 3 

DISCUSSION OF RESULTS 


In this section, first results are presented for the case without jet, followed by the results 
for the case with jets. As mentioned earlier, complete analysis is carried out using a well known 
and widely used commercial CFD code, Fluent. In CFD analysis, the need for validating the 
results by comparing with experimental results or by other means is well established. Many times, 
this process becomes a very difficult task due to lack of experimental and numerical results. In the 
present study, results are compared with existing data or with the data obtained using similar 
physical configurations, but not exactly the same. In the latter case, effort is directed to justify 
when such comparisons are made. 

It is important to note some typical observations regarding flow past circular and rectangu- 
lar cylinders based on the information available in the literature. One important similarity between 
these two configurations is the vortex shedding phenomena observed in both cases. The important 
difference between the two configurations is the point of flow separation. In case of rectangular 
cylinders, the flow separates at the leading edge, even for very low Reynolds number [5], The 
point of flow separation, in case of circular cylinders, is highly dependent on the free stream Rey- 
nolds number. The point of separation moves away from the wake and towards the upstream with 
increase in Reynolds number [14], These well established observations can be utilized for validat- 
ing numerical results if necessary. 

The Strohaul number, base pressure coefficient, and lift and drag coefficients, are com- 
puted for the case with no jets for different free stream Reynolds numbers. Few cases are also con- 
sidered using different aspect ratios. In the case of situation with jets, drag coefficients are 
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evaluated for different free stream Reynolds numbers. The effect of introducing a jet on drag is 
analyzed. 

3.1 Case 1: No Jet 

Unsteady state computation is carried out when no jet is considered at the rear stagnation 
region. As mentioned in the previous section, the values of time steps for different inlet conditions 
(Reynolds numbers) are chosen after several trial tests. The time steps chosen for different Rey- 
nolds numbers is listed in the Table 3.1. Free stream Reynolds number is varied from 1 X 10 4 

to 1.25 X 10 5 . Size of the time step decreases with the increase in free stream Reynolds number. 
The physical reason for this trend is, due to the increase in Reynolds number, the change in flow 
takes place rapidly. Consequently, a small time step is needed to capture the physical processes 
more accurately. 

The Strouhal number variation with the free stream Reynolds number is listed in the Table 
3.2. It is a general practice to numerically compute the Strouhal number based on the fluctuations 
of pressure or lift coefficient [7]. Here, the Strouhal number is estimated based on the time of one 
cycle of the pressure fluctuation. Figures 3. 1-3.5 show the pressure fluctuation with time for dif- 
ferent Reynolds numbers considered in this study. The results were computed on the center point 
of top edge of the rectangle with aspect ratio equal to unity. From these plots, frequency needed to 
estimate the Strouhal number is obtained. The results for pressure fluctuation with time for differ- 
ent Reynolds numbers are shown in Fig. 3.6. The frequency of pressure fluctuation increases with 
the increase in the free stream Reynolds number. The variation of the Strouhal number with the 
free stream Reynolds number is shown in Fig. 3.7. For the range of the free stream Reynold num- 
ber used in this study, the value of the Strouhal number remained constant at about 0.12. Atsuchi 
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Table 3. 1 Size of time step for different free stream Re numbers 


Reynolds 

number 

1 X 10 4 

2.5 X 10 4 

50 X 10 4 

75 X 10 4 

1.25 X 10 5 

Time step 
size, in Sec. 

0.1 

0.05 

0.02 

0.02 

0.01 
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Table 3.2 Variation of Strouhal number with free stream Re number( AR = 1 ) 


Reynolds 

number 

1 X 10 4 

2.5 X 10 4 

5.0 X 10 4 

7.5X10 4 

12.5 X 10 5 

Strouhal 

number 

0.1187 

0.1194 

0.1179 

0.1114 

0.1037 
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Pressure, N/m 



Fig. 3.2 Time history of pressure (Re = 2.5 X 10 4 ) 
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- 74.0 



Time, Sec 

Fig. 3.5 Time history of pressure (Re = 1.25 X 10 5 ) 
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Pressure, N/m 



Time, Sec 


Fig. 3.6 Time history of pressure (Re = 1 X10 4 - 1.25 X10 5 ) 



Strouhal number 



Fig. 3.7 Variation of Strouhal number with Reynolds number 




[5] measured experimentally the Strouhal number in the range of the free stream Reynolds num- 
ber between 70 and 2 X 10 4 . In this range of Reynolds numbers, it was reported that the Strouhal 
number of the square cylinder showed slight and continuous change around a constant value of 
0. 125. The present result also compares well with the measured value of the Strouhal number by 
Bearman and Truman [3]. Once again, the comparison is restricted to the lower Reynolds number 
range due to non-availability of data at higher Reynolds number range. 

The Strouhal numbers are also estimated for different aspect ratios. Figure 3.8 shows the 
fluctuations of pressure for aspect ratio equal to 0.8 and 1.5. Table 3.3 presents the estimated val- 
ues of the Strouhal number for different aspect ratios when Reynolds number is equal to 1 X 10 4 . 
The value of the Strohaul number nearly remained constant with the aspect ratio considered in 
this study. 

The base pressure coefficient along the down stream face (base) is estimated next. Here, 
the base pressure coefficient is obtained by taking average over a number of cycles, once system- 
atic periodic oscillation is established. The variation of the base pressure coefficient with Rey- 
nolds number is shown in Fig. 3.9. Experimental value obtained by Bearman and Truman [3] 
matches well with the estimated value of the base pressure coefficient. Once again, comparison is 
restricted to the lower Reynolds number range. 

Next, drag and lift coefficients are plotted versus time for different Reynolds numbers con- 
sidered in this study. These plots are shown in Figs. 3.10-3.19. The systematic oscillations can be 
seen in these plots. Average lift and drag coefficients are listed in Table 3.4 for different Reynolds 
numbers. Comparison of estimated average drag coefficients is also made with experimentally 
measured values at the lower Reynolds number range which is shown in Table 3.5. The value of 
drag coefficient given by Nakaguchi et al. is taken from the paper by Beraman and Truman [3] 
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Pressure, N/m 



Fig. 3.8 Time history of pressure for different aspect ratios (Re = 1 X10 4 ) 
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Table 3.3. Variation of Strouhal number with aspect ratio ( Re = 1X10 4 ) 


Aspect ratio 

0.8 

1.0 

1.5 

Strouhal 

number 

0.1311 

0.1187 

0.1338 
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Fig. 3.9 Variation of base pressure coefficient with Reynolds number 



Drag coefficient 


AR= 1.0 



Fig. 3.10 Time history of drag coefficient (Re = 1 X 10 4 ) 
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Drag coefficient 



Fig. 3.13 Time history of drag coefficient (Re = 75 X 10 4 ) 
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Lift coefficient 



Fig. 3.16 Time history of lift coefficient (Re = 2.5 X 10 4 ) 
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Lift coefficient 




Lift coefficient 
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Lift coefficient 



Fig. 3. 19 Time history of lift coefficient (Re = 1.25 X 10 5 ) 
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Table 3.4 Variation of C d and Q with Re number 


Reynolds 

number 

1 X 10 4 

2.5 X 10 4 



1.25 X 10 5 

Drag coefficient 

1.825 

1.832 

1.855 

1.834 

1.50 

Lift coefficient 

0.162 

0.123 

0.145 

0.1174 

0.118 
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Table 3.5. Comparison of C d with experimental value( Re = 1 X 10 4 ) 



Present study 

Beraman and 
Truman [3] 

Nakaguchi et al. 

Drag coefficient 

1.825 

2.00 

1.90 
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Finally, streamlines are plotted at the maximum, middle, and minimum lift coefficient 
points after periodic oscillation is established. Figure 3.20 shows the streamline pattern for Rey- 
nolds number of 1 X 10 4 . Flow separates immediately at the leading edge, even for the lowest 
Reynolds number considered in this study. Similar results are presented in Figs. 3.21-3.23 for 
other Reynolds numbers. As can be seen from these plots, the streamline pattern changes with 
time. The streamline pattern also changes with increase in Reynolds number. 

3.2 Case 2 : With Jet 

The important reason for a jet in the rear stagnation region is for reducing the drag. A 
recent experimental and numerical study by Atsuchi et al. [10] gives information on various 
aspect of fluid dynamics, in the case of a circular cylinder. The cross flow over a two-dimensional 

circular cylinder at Re = IX 10 5 was used as a “bench mark'’ test. Since, no information is avail- 
able in case of rectangular cylinders, it is reasonable to make some judgements based on this work 
as the overall aspects of fluid dynamics is very much similar. Only effect of free stream Reynolds 
number is investigated. Blockage ratio, aspect ratio as well as jet blowing momentum coefficient 
C m are unchanged in this study. 

In case of circular cylinders, when the jet blowing momentum coefficient is lower than 
0.0005, no change is seen in the drag coefficient. Drag decreases for C m between 0.001 and 0.02 
[10]. It was also noted that further increase in C m can actually increase the drag coefficient. 
Although, it is not possible to draw the same conclusion with complete confidence, in case of rect- 
angular cylinders, it was decided to use intuitively a value of 0.01 for C m . As pointed in this sec- 
tion, this value of C m helped in drag reduction. The study is conducted using the free stream 

Reynolds number ranging from 1 X 10 4 to 1.25 X 10 5 . Also in this case (with jet), steady state 
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assumption is made. Drag coefficients are listed in the Table 3.6 for different Reynolds numbers 
for both cases, i.e, with jets as well as without jets. It is clearly seen that the drag is reduced when 

jet is introduced at rear stagnation region. Only in the case of Re = 1.25 X 10 5 , the change in drag 
is negligible. 

Figure 3.24 shows streamlines plots for different free stream Reynolds numbers with C m = 

0.0 1. The streamline pattern becomes more symmetric when jet is introduced in the rear stagna- 
tion region which could be the reason for reduction in drag coefficient. 

Further, numerical experimentation is carried out to investigate the effect of increasing C m 
on drag coefficient. Drag coefficients are listed in the Table 3.7 for different values of Cm when 
Re = 1.25 X 10 5 . The value of drag coefficient is not affected when C m =0.0005. It decreases 
when the value of C m is increased to 0.01, but further increase in C m actually increases drag coef- 
ficient. Figure 3.25 shows streamlines plots for the different values of C m at Re = 1.25 X 10 5 . The 
position and size of recirculation zone clearly changes with the value of C m which is the main rea- 
son for affecting the drag coefficient. 
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Table 3.6 Value of Cj with jet and no jet for different free stream Re number 


Reynolds 

number 

1 X 10 4 

25 X 10 4 

50 X 10 4 

75 X 10 4 

1.25 X 5 

Drag coefficient 
(without jet) 

1.825 

1.832 

1.855 

1.834 

1.50 

Drag coefficient 
(with jet, C m =0.01) 

1.557 

1.596 


1.534 

1.55 
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(c) AtRe= 1.25 X 10 5 

Fig. 3.24 Streamlines from steady computation ( C m = 0.01 ) 
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Table 3.7 Value of C d for different C m at Re = 1.25 X 10 5 


C m 

0.0005 

0.01 




1.83 

1.55 

1.89 

1.90 
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Chapter 4 


CONCLUDING REMARKS 

Numerical investigation of flow across rectangular cylinders is conducted with and with- 
out jets in the range of Reynolds numbers between 1 X 10 4 to 1.25 X 10 5 . The Strouhal number, 
average base pressure coefficient, lift and drag coefficients are numerically estimated in the 
absence of jet. When jet is considered, possibility of drag reduction is examined. Drag coefficient 
decreases initially with an increase in jet momentum coefficient C m but beyond a value of 0.01, 

C d actually inceases. 
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